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Abstract. The paper is devoted to the understanding of the infrared emission of nuclear regions in galaxies, a) ISO 
data of NGC 1808 are presented: spectro-photometry from 5.1 to 16.4 /im; a 25" x 25" map at 6[im and 170/^m 
photometry, b) The data are complemented by a polarization measurement at 170/im (2.5±0.4% at position angle 
94° ±5°) and a map at 6/im. In the map, the degree of polarisation goes up to 20% in the outer regions. We argue 
that the polarisation is produced by emission of big grains and exclude very small grains and PAHs or scattering 
and extinction, c) The mid infrared spectrum shows, beside the main emission bands, a so far unknown plateau 
of PAH features in the > 13 /j,m region, d) The total spectrum can be fit under the assumption of optically thin 
emission. However, such a model fails to reproduce the 25^m point and implies that the mid infrared is due to 
very small grains and PAHs. These particles would then also have to be responsible for the 6/im polarisation, 
which is unlikely, e) To avoid these difficulties, we successfully turn to a radiative transfer model whose major 
feature is the existence of hot spots produced by the dust clouds around OB stars. We demonstrate the decisive 
influence on the mid infrared spectrum of both the PAHs and the hot spots. 
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, ^ 1. Introduction 



The nucleus of a galaxy is usually deeply embedded in 
dust and can therefore only be probed by infrared (IR) 
observations. In this paper, we present an infrared polar- 
isation and a spectro-photometric study of the starburst 
galaxy NGC1808. The new data are explained by a con- 
sistent model for grain heating and cooling. We include a 
discussion of the various infrared emission bands (IEB). 
Understanding the mid IR emission is of particular inter- 
est since the ratio of mid IR over total luminosity is a 
diagnostic for estimating the relative contribution of star- 
burst versus AGN activity (Rigopoulou et al. 1999, Clavel 
et al. 2000). 

Another important and so far unanswered question is 
how much of the galactic infrared emission is polarised. 
The observed emission in galaxies is always averaged 
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over a large volume and, in case polarisation is detected, 
should be linked to large scale alignment of the dust. 
Polarised emission from interstellar grains can only be 
detected at infrared and submillimeter wavelengths and 
has been observed in various galactic sources, such as 
the Galactic Centre (Hildebrand 1989), molecular clouds 
(Novak et al. 1989, Clemens et al. 1999), star forming re- 
gions (Gonatas et al. 1990, Leach et al. 1991, Smith et 
al. 2000), protostars (Siebenmorgen & Kriigel 2000) and 
prestellar cores (Minchin et al. 1995, Ward-Thompson et 
al. 2000), and in the nucleus of the starburst galaxy M82 
(Greaves et al. 2000). 

The galaxy NGC1808 is at a distance of 11.1 Mpc 
(Reif et al. 1982), where 1" corresponds to 53 pc. It is 
classified as type Sbc (Sandage & Tammann 1987) and 
seen fairly edge-on. The moderate far IR luminosity of a 
few 10 10 L Q is ascribed by Krabbe et al. (1994) to young 
stars formed in recent star formation. For the central 
1.5 kpc, Koribalski et al. (1993) suggest a total gas mass 
(HI plus molecular gas) of 4.2xl0 9 M Q , and Dahlem et 
al. (1990) propose 2xl0 9 M© of molecular gas. In the op- 
tical, NGC1808 reveals a disk with outer spiral arms. Dust 
filaments are seen in the North-East and seem to emerge 
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from the nucleus (Veron-Cetty & Veron 1985). There is 
evidence of a central bar of 500 pc length, an inner ring- 
like structure (Tacconi-Garman et al. 1996), and a number 
of compact radio (Saikia et al. 1990) as well as near IR 
sources (Kotilainen et al. 1996). 

In Sect. 2, we present the observations followed by a 
discussion of the infrared polarisation (Sect. 3). A model 
for astronomical PAHs is outlined in Sect. 4. The infrared 
spectral energy distribution is first discussed assuming op- 
tically thin emission (Sect. 5); the model is refined in 
Sect. 6 to include proper radiative transfer. Results are 
summarised in Sect. 7. 

2. Observations 

ISOCAM (Cesarsky et al. 1996) polarimetric imaging on 
NGC1808 was set up according to the observing template 
CAM05 (Siebenmorgen et al. 1996). The corresponding 
ISO observation numbers, TDT, are: 87501602, 87501703, 
87501804, 87501905, 87502006, 87502107 and 87502208. 
We used the 6/im filter lw4 and a 3 x 3 raster with a raster 
step size of 39" and the 3" lens. The read out time of each 
exposure was 2.1s. After a background measurement, we 
took 15 exposures on the source to stabilize the detec- 
tor and then a raster through the free hole of the en- 
trance wheel. At each raster position 15 exposures were 
read out. This procedure was carried out for each of the 
three polarisers. The polariser rasters were repeated in 2 
cycles. The data were reduced with the ISOCAM inter- 
active analysis system (CIA version 4.0, Ott et al. 1996). 
Only basic reduction steps, such as dark current subtrac- 
tion, removal of cosmic ray hits or transient correction, are 
applied to the data. They are described together with cal- 
ibration uncertainties in Siebenmorgen et al. (2000). The 
coadded images at each raster position were projected on 
the sky and corrected for field distortion to derive the fi- 
nal mosaics. The mosaics have a 168" x 168" total field of 
view and a pixel scale of 1". The polarised signal is found 
to be consistent between the cycles, and the average of all 
cycles gives the final mosaic image for each polariser. Our 
polarisation image is shown in Fig. |l| after correction for 
instrumental polarisation (Siebenmorgen 1999). 

From the ISO archive we retrieved spectro- 
photometric-imaging data taken with the circular 
variable filter (CVF) of ISOCAM (TDT: 69702381). 
For each CVF step between 16.33/zm down to 5.14/zm, 
11 exposures of 2.1s integration time using the 3" lens 
were read out (Siebenmorgen et al. 2000, for more infor- 
mation on ISOCAM data handling). After background 
subtraction, we simulate multi-aperture photometry on 
the images. The spectrum, as derived for a 25" aperture 
centred on the brightest pixel, is shown in Fig. H and 
Fig. ||. One notices strong IEBs centered at 6.2, 7.7, 8.6, 
11.3 and 12.8 ^m. There is a weaker but definite broad 
band emission feature centered at 7.0/xm. A detection of 
the 7.0/im band is reported for HD9730 by Siebenmorgen 
et al. (1998) and in other galaxies by Helou et al. (2000). 
In addition there is a plateau of fainter bands beyond 




x offset (arcsec) 

Fig. 1. Polarisation image of NGC1808 at 6/im (ISOCAM 
filter lw4) . The observed E-vectors and the contours of the 
total emission are shown in detector coordinates. The ce- 
lestial North-East orientation is indicated. Contours range 
from 1.28 to 8.21 in steps of 0.77 mJy/arcsec 2 . The effec- 
tive resolution is 5". 

13/.tm and most likely we see the rise towards the 16.4 
//m band (Moutou et al. 2000), which is just outside the 
scanned wavelength region. Below we combine the main 
bands as well as the fainter bands and newly discovered 
IEBs with measured laboratory modes of polycyclic 
aromatic hydrocarbons (PAH) and present quantitative 
model fits. 

Linear polarisation at 170 /im was detected with the 
ISOPHT C200 2x2 array (Lemke et al. 1996) by mea- 
suring several cycles with the 0°, 120° and 240° polarisers 
(TDT: 85601702, 85601803, 85601904 and 85602005). The 
cycles were followed by an open (sky) measurement and 
one of the internal reference calibrators (Klaas et al. 1999). 
We commanded 4 pointings such that each of the 4 detec- 
tor pixels measured the polarisation of the centre position 
of the galaxy in succession. In addition, the observations 
on the source were bracketed by two background obser- 
vations to determine the zero signal level and to monitor 
instrumental variations. The sky measurements were pro- 
cessed and calibrated using the standard procedures for 
photometry (Laureijs et al. 2000) with the ISOPHT inter- 
active analysis software (PI A version 8.1). After coaddi- 
tion of the 4 pointings of the C200 array, we obtained a 
small map which indicated that the brightness distribution 
of the galaxy is not significantly different from that of a 
point source. We conclude that NGC1808 is not resolved 
with ISOPHT at 170 fiva. The galaxy is much brighter 
than the background and has a flux i^iro^m = 104±10 Jy, 
most of which comes from the central pixel. The estimated 
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error of 10% is based on the ISOPHT absolute photomet- 
ric calibration uncertainty. 

Since the galaxy is a point source at 170/im, we de- 
termined the polarisation only for the four pixel-position 
combinations centred on the galaxy. The polarisation mea- 
surements of the other positions in the map were not used 
as they were pointing on the wings of the beam profile of 
the galaxy. Due to the uncertainties in the beam profiles 
of each detector pixel, we are not sure whether our values 
for instrumental polarisation apply to measurements for 
which the detector pixel is not centred on the (polarised) 
point source. The processing of the signals were carried 
out according to the procedure described in Laureijs & 
Klaas (1999); we also used their instrumental polarisation 
values. The resulting fractional polarisation degree (p) and 
polarisation angle (8) with respect to the equatorial North 
per C200 detector pixel are presented in Table g- The po- 
larisation vector obtained from pixel 3 is significantly dif- 
ferent from the other pixels and has a somewhat lower 
accuracy. The last line in Table @ gives the polarisation 
derived from the average of the Stokes vectors of all pixels 
after subtraction of the instrumental polarisation. We did 
not include possible systematic uncertainties introduced 
by the instrumental polarisation of the pixels. We estimate 
that they amount to less than 1% in the final polarisation 
degree. 

Table 1. Polarisation of the centre position of NGC1808 
at 170/im for each pixel of the ISOPHT C200 array. Due to 
the 180° ambiguity, the angles in the first two quadrants 
were chosen. 



pixel 
number 


P 

[%] 


9 
[degrees] 


1 
2 
3 

4 


4.1 ± 0.9 
5.3 ± 1.2 

2.2 ± 1.1 
5.0 ± 1.0 


75 ± 6 
116 ± 6 
3 ± 14 
114 ± 6 


mean 


2.5±0.4 


94±5 



3. First estimates 

As the source is not resolved by ISO at 170/im, its 
far IR size must be smaller than the diameter of the 
first Airy ring (140"). Scans at lOOfim. from the KAO 
(Smith & Harvey, 1996) indicate a maximum source size 
of ©source = 36" compatible with the maps at 6/xm (this 
paper), and at 8.4 GHz and in Br7 (Kotilainen et al. 1996). 
To get an estimate of the average optical depth, we put 
F\ 7o M m = ^ • Bir 0fim (T) • Ti7o Mm and assume homoge- 
neously distributed dust within il = 2.39 • 10~ 8 steradian 
and at temperature of 33K, as is typical for starburst 
galaxies (Siebenmorgen et al. 1999). In this way, one gets 
Ti70/im ^ 6.5 ■ 10~ 3 in rough agreement with the estimate 
by Smith & Harvey (1996). By applying a standard dust 
model (Krugel & Siebenmorgen 1994b), this translates 



into a visual optical depth ry ~ 2000 • r^o^m < 13 or 
r 6Mm ~ r v /60 < 0.2. 

The fractional polarisation pi7o Mm = 2. 5±0. 4% implies 
a total polarised flux of ~ 2.5 Jy. The extrapolation of the 
(non-thermal) 20 cm radio flux of 520 ± lOmJy with a 
^~°- 7 -law to 170/im gives only 4mJy. Therefore the far IR 
polarisation of NGC1808 cannot contain any synchrotron 
contribution, but must be entirely due to thermal emission 
from dust grains aligned by a magnetic field. 

The polarisation angle #i7o M m = 94° ± 5° indicates 
the direction of the E-vector, which is perpendicular to 
the interstellar magnetic field. This result differs from 
the strength and field direction derived from the non- 
thermal radio emission. Averaging the central 90" of the 
6 cm map by Dahlem et al. (1990), the angles differ by al- 
most 70°. However, both observations may probe regions 
in the galaxy at different vertical heights where the mag- 
netic field has a different horizontal component. A good 
example for such a case is presented by M51 (Fig.l of 
Berkhuijsen et al. 1996). 

Our 6/mi ISOCAM map of NGC1808 also depicts the 
central region. Its size down to the lowest contours equals 
23" x 15" and the integrated flux is 1130±60mJy. A com- 
parison with the Br 7 emission (Kotilainen et al. 1996), 
which follows the OB stars, shows that we also trace the 
starburst at 6/^m. The starburst occurs mainly along the 
major axis of Fig. El Polarisation was also detected in the 
ISOCAM map. The polarisation is strongest (up to 20%) 
in the outer regions and weak in the centre of the star- 
burst. Along the major axis and the central region of the 
galaxy where the Br7 emission and the near IR and radio 
hot spots are concentrated, the degree of polarisation is 
quite small. Because the E-vectors have a uniform direc- 
tion over distances of 10", which is greater than the instru- 
mental angular resolution, the field itself must be ordered 
over such a scale, corresponding to 500pc or more. 

Extrapolating the 20 cm total radio flux towards the 
mid IR, we can exclude any net polarisation from non- 
thermal emission. When we integrate at 6/xm over the 
whole nucleus, we do not detect any net polarisation, but 
obtain only a 3a upper limit of 5%. 

In the optical (V band), the degree of polarisation 
amounts to a few percent (Scarrot et al. 1993), which is 
much smaller than in the mid IR. The polarisation vectors 
at V are quite regular and lie perpendicular to the radial 
vectors towards the center of the galaxy. In the inner part 
(750pc) of the galaxy, the polarisation is reminiscent of a 
reflection pattern powered by the nucleus. 

The mid IR polarisation, on the other hand, is also 
fairly regular, but neither circular nor centre-symmetric. 
This argues against scattered light from the galactic nu- 
cleus. Moreover, the grains are tiny compared to the wave- 
length so that their scattering cross section must be very 
small. Therefore, the mid IR polarisation mechanism is 
different from the optical and is either caused by extinc- 
tion or emission of aligned dust. In any case, some of the 
interstellar grains must be elongated and aligned, a re- 
quirement not necessary for polarisation by scattering. 
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Table 2. Astronomical PAH 



(1) 


(2) 


(3) 


(4) 


(5) 


Wavelength 


Damping Constant 


Integrated Cross Section 


Mode 


Remarks 


Ao 


7 


0"int 






/im 


io 1 ^- 1 


10~ 22 cm 2 /im 






3.3 


20 


12 


C-H stretch 


no data 


5.2 


12 


1.1 


C-C vibration 


tentative 


6.2 


14 


21 


C-C vibration 




7.0 


5.9 


12.5 


C-H? 




7.7 


22 


55 


C-C vibration 




8.6 


6 


35 


C-H in-plane bend 




11.3 


4 


36 


C-H solo out-of-plane bend 




11.9 


7 


12 


C-H duo out-of-plane bend 




12.8 


3.5 


28 


C-H trio out-of-plane bend 




13.6 


4 


3.7 


C-H quatro out-of-plane bend 


tentative detection 


14.3 


5 


0.9 


C-C skeleton vibration 


first extragalactic detection 


15.1 


3 


0.3 


C-C skeleton vibration 


tentative detection 


15.7 


2 


0.3 


C-C skeleton vibration 


first extragalactic detection 


16.4 


3 


0.5 


C-C skeleton vibration 


data do not cover full band 


18.2 


3 


1.0 


C-C skeleton vibration 


no data 


21.1 


3 


2.0 


C-C skeleton vibration 


no data 


23.1 


3 


2.0 


C-C skeleton vibration 


no data 



(1) Center wavelength Ao. 

(3) Cross sections intregrated over the band. They are determined from our radiative transfer model and given per H atom for 
C-H vibrational modes, and per C atom for C-C modes. 

(4) In the assignment of the weak bands, we follow Moutou et al. (1996) and Pauzat et al. (1997). 



If the 6/im polarisation at a 10% level or more is due to 
extinction, the associated optical depth at 6/^m must be 
greater than unity (Siebenmorgen & Kriigel 2000, Fig. 11) 
and four or more at 10/Ltm. This would imply a very strong 
silicate absorption feature which is not seen. We therefore 
favor the idea that the polarisation is due to grain emis- 
sion, as in the far IR. Note that the -B-vectors of the 
aligning field for the two scenarios are perpendicular to 
each other. 

The polarised mid IR emission of the dust must have one 
of the following origins: 

1. It may come from PAHs heated by the interstellar radia- 
tion field (ISRF) in the galactic nucleus (which is certainly 
much stronger than the ISRF near the sun). Indeed, our 
ISOCAM mid IR spectrum (Figs. 4 and 5) shows strong 
IEBs. The orientation of a PAH along an interstellar mag- 
netic field works only if it has a magnetic moment and if 
the timcscalc for alignment is smaller than for disorder. 
Although carbon and hydrogen atoms are not paramag- 
netic, Rouan et al. (1992) propose a number of alterna- 
tive ways how to endow PAHs with a magnetic moment. 
However, their conclusion is that PAHs, whether ionized 
or not, cannot align. Therefore we do not favour PAHs to 
produce the observed polarisation. 



2. Small silicates are paramagnetic and consequently could 
align. But to account for the polarisation they must also 
contribute sufficiently to the emitted flux. Calculation 
of the emission from various types of transiently heated 
grains in reflection nebulae (Siebenmorgen et al. 1992, 
1998; Li & Draine 2001b) show that small silicates are 
not efficient radiators at 6/im. 

3. Large grains with sizes a > 100 A, must have temper- 
atures above ^200 K to emit substantially at 6/xm. The 
heating by the ISRF is insufficient and the only possibil- 
ity for them to get so hot is close to OB stars (in so called 
hot spots, see Fig. |fj and Sect. |fj). Large grains are our 
favored alternative. 



4. The dust model 

4.1. Large and small grains 

The dust model adopted in the following two sections con- 
sists of amorphous carbon particles (abundance in solid 
[C]/[H] = 2.9 • 10- 4 ) and silicates ([Si]/[H] = 3 • 10" 5 ) 
with bulk density 2.5g/cm 3 . Their radii range from 100 to 
2400 A with an a -35 size distribution (Mathis et al. 1977). 
Absorption and scattering cross sections are calculated 
from Mie theory employing optical constants from Zubko 
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Fig. 2. Infrared absorption coefficient k\ of astronomical PAH. There are two populations, one with iVc = 30 C 
atoms, the other with Ac = 500. Both populations have an hydrogenation parameter /h/c— 0-2 and an abundance 
such that for 10 5 protons in the gas phase one finds one carbon atom in each PAH population (see also Table ||). 



et al. (1996) for amorphous carbon and Laor & Draine 
(1993) for silicates. 

There is another dust population of very small 
graphites and very small silicates with radii 5 A< a <80 A 
and size distribution n(a) oc a~ 4 whose temperatures fluc- 
tuate. Cross sections are again from Mie theory optical 
constants for graphites come from Draine & Lee (1984). 
The abundances of the small graphite and small silicate 
grains are typically 5-10% that of the large carbon and 
large silicate grains. The assumed abundances imply a 
gas-to-dust mass ratio of ^140. 



4.2. Astronomical PAH 

For simplicity, we consider only two kinds of PAHs: small 
ones with Nq = 30 carbon atoms and big ones with Ac = 
500. Both have a hydrogenation parameter, which is the 
ratio of hydrogenation to carbon atoms, of /h/c = 0-2 
and a carbon abundance of Y"<? AH = 3% relative to the big 
grains. We calculate absorption cross sections in the UV 
and visible following Schutte et al. (1993) and disregard 
the differences between neutral ond ionized PAH species 
(De Frees et al. 1993, Onaka et al. 1999, 2000, Chan et 
al. 2000, Uchida et al. 2000). For the infrared emission 



bands of the C-H and C-C vibrational modes we assume 
Lorentzian profiles: 



N-a h 



"{V 



C TT 2 (v 2 -4) 2 + ( 1 v/2) 



(1) 



where 7 is the damping constant, a lrA the integrated 
cross section and N the number of carbon or hydrogen 
atoms of the PAH in the particular band at center fre- 
quency vq. Model parameters of the PAH resonances are 
listed in Table ||. Altogether, we consider 17 PAH bands 
most of which are actually detected in NGC1808, except 
the one at 3.3/im and those whose center wavelength is 
greater than 16.3/im. 

The assignment of the feature at 7/tm is unclear, it 
could either be a C-C stretch or due to C-H (Allamandola 
et al. 1989). There may also be some blending with gas 
emission lines (e.g. [Ar II] line at 6.983 /im) which can- 
not be resolved with the CVF. The bands with 13.6 /im 
< Ao < 16.3/im are seen for the first time in an extra- 
galactic object, but are known from Milky Way sources 
(Van Kerckhoven et al. 2000, Moutou et al. 2000, Hony et 
al. 2001). 

In the far IR, properties of PAHs are not well con- 
strained and we assume the continuum to fall off propor- 
tional to A -12 for A > 24/im. Our infrared absorption 
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Fig. 3. Absorption coefficient n\ from the UV to the near IR for the dust model of Section 4.2: Sum of large carbon and 
large silicate grains (dash-dot), small graphites (dashed) and small silicates (full line). PAHs of 30 C atoms (dotted) 
with a cutoff at ~ 0.26/im. The absorption coefficient of the PAHs with Ac = 500 is the same, except that their cutoff 
is shifted to 1/mi. 



cross sections are displayed in Fig. g. They lie within the 
range of lab measurements (Moutou et al. 1996, Pauzat et 
al. 1997, Schutte et al. 1993) and they may be compared to 
other models, for example, by Li & Draine (2001). Fig. || 
shows the contribution of the dust populations to the total 
absorption cross. 

5. An optically thin model 

As a first attempt, we model the spectral energy distribu- 
tion of NGC1808 under the assumption that the emission 
is optically thin and the dust mixture is exposed to diluted 
black-body radiation of mean intensity 



J v 



B„(T) 



(2) 



where the dilution factor, a, is a free parameter. The 
temperature T reflects the effective temperature of the 
stellar population in the nucleus; we use T = 20000 K. In 
the model depicted in Fig. |], a — 2.2-10~ 15 . In this model, 
we have slightly modified a few integrated cross sections 
of the PAHs with respect to the numbers of column (3) 



in Table @, but never more than 30%. The total gas mass, 
employing the abundances of the dust components as de- 
scribed above, equals Af gas = 4 • 1O 9 M . 

Overall the fit to the data in Fig. |j is good and for the 
PAH bands even excellent. A deficit of the model around 
5.7/im could possibly be due to yet unresolved band emis- 
sion. Such a band is apparent in other spectra, for exam- 
ple, of the Orion bar and BD+30 3639 (Fig. 10 by Schutte 
et al. 1993). The 5.2/zm is not well constrained by the data. 
The 7.0/im IEB is detected. We point out the blending of 
the 12.0 /im PAH band with the adjacent main features. 
In the optically thin model, the 6/xm radiation, which is 
further investigated in our polarisation study, is due to 
PAHs. 



6. Radiative transfer: Hot spots 

In the previous section we have matched the infrared spec- 
tral energy distribution of the nucleus of NGC1808 by 
optically thin dust emission and reasonable assumptions 
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Fig. 4. Top: Model spectrum of NGC1808 for optically thin emission. The large grains emit only in the far IR (dash- 
dot line), small graphites (dotted) and small silicates (dashed) do not contribute significantly to the total emission 
(full line). Data points from IRAS, ISOCAM, ISOPHT and Glass (1976). In the mid IR, the emission is dominated 
by PAHs. Bottom: Comparison of optically thin model after convolution with the ISOCAM CVF spectral response 
function (full line) with the ISOCAM spectrum. 
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Fig. 5. Top: Radiative transfer model spectrum of the nucleus of NGC1808. We show the model flux as received in a 
25" (dots) and a 100" aperture (solid). Data points arc shown for ISOPHT (100" resolution), IRAS, ISOCAM (25" 
aperture, barely visible), and in the near IR (25" aperture) by Glass (1976). For model description see text. Bottom: 
Blowup of the near and mid IR data and the 25" aperture model. 
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Table 3. Parameters of the radiative transfer model. 



Parameter 


NGC 1808 


Distance 


11.1 Mpc 


Extinction 




Ay to center (r=0) 


4.3 mag 


Luminosity 




of all giants Lg 


5.5 • 10 10 L 


of all OB-stars Lob 


0.5 ■ 10 10 L 


of one OB star Lhs 


10 4 L Q 


Temperature 




of giants 


5000 K 


of OB stars 


30000 K 


Outer radius 




of nucleus Rn 


3000 pc 


of giants 


3000 pc 


of OB stars Rob 


350 pc 


Density 




of nucleus 


r < Rob'- const 




-Rob < r < _Rn : oc r _1 


of giants 


OCT-" 1 ' 8 


of OB stars 


oc r _1 


in hot spots 


10 4 H atoms per cm 3 


Small PAHs 




abundance Yq h 


6% 


size Ac 


50 
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about the exciting radiation spectrum. Such a model gives 
a good match to the CVF and far IR data, however, 

1. it cannot account for the 25^m emission; 

2. it gives no information of the structure within the nu- 
cleus nor on the observed IR sizes; 

3. it incorrectly assumes zero optical depth within the 
nucleus; 

4. its 6/im emission is due to PAHs. As PAHs are unlikely 
to be good polarisers, this is in contradiction to our IR 
polarimetric results. 

We therefore include the effects of radiative transfer 
in the nucleus of the galaxy. Such computations have 
been carried out in various approximations by a num- 
ber of authors (Rowan- Robinson & Crawford 1989, Pier 
& Krolik 1993, Granato & Danese 1994, Efstathiou & 
Siebenmorgen 1995, Silva et al. 1998, Siebenmorgen et 



al. 1999b, Efstathiou et al. 2000, Ruiz et al. 2001). 
The procedure applied here is described by Kriigel & 
Siebenmorgen (1994). The model assumes spherical sym- 
metry which seems justified as the mid IR continuum im- 
age of NGC1808 (Fig. pi) indicates an oblate structure with 
ratio of major-to-minor axis of less than 3 : 2. It has the 
following characteristics: 

— There are two stellar populations both distributed 
within the nucleus: a) giants of 5000K surface temper- 
ature with a radial density distribution that falls off 
from the center to the edge of the nucleus like r^ 18 . 
b) OB stars (T OB = 30000K, L OB = 10 4 L Q ) each of 
which is surrounded by a spherical dust cloud. Their 
distribution drops inversely proportional to the radial 
distance from the nucleus. 

- There are hot spots (Kriigel & Tutokov 1978). The 
heating of the dust in the galactic nucleus does not 
only change on large scales, but also locally. A grain 
at a galactic distance R, as far as possible away from 
any star, gets less heated than one at same galactic 
radius R, but close to a star. The hot spots are formed 
by the dusty envelopes around the OB stars. 

In computing the radiative transfer in the hot spots, we 
assume that they are spheres. On the surface of the spheres 
the dust is heated equally by the ISRF of the galactic 
nucleus and by the central OB star. Consequently, the 
size of the hot spots is smaller in the center and increases 
in the outer regions of the nucleus. 

Note that the emission of the hot spots is not just 
added to the emission from the rest of the nucleus, 
but is included in a self-consistent way (see Kriigel & 
Siebenmorgen 1994). The main effect of the hot spots is 
to drastically enhance the mid IR flux. 

The model nucleus of NGC 1808 has a total luminosity 
Ltot = 5.6 • 10 10 L o , a gas mass M gas = 1.2 • 10 9 M Q , not 
terribly different from CO estimates (Dahlem et al. 1990), 
and a radius R ~ 3kpc. The dust density is constant in the 
inner parts of the nucleus where the OB stars are located 
and declines smoothly like p(r) ex r _1 further out. For the 
distribution of the OB stars, we take the FWHM of the 
6/im image as a guide and adopt an outer radius Rob = 
350 pc. From fitting the near IR data, we hnd that the near 
IR flux is mostly due to giants and put Lq = 5.5 • 10 10 L . 
The OB stars have a total luminosity Lob = 0.5 • 10 10 L@, 
each individual OB star or hot spot has Lhs = 10 4 L Q . 
The density in the hot spots is constant, nns(H+H2) = 
10 4 H atoms per cm 3 . Model parameters are summarised 
in Table |. 

Results are shown in Fig. g. Because of the hot spots, 
the 25/im IRAS point is matched, which was not possible 
in the optically thin model of Fig. [|, and the agreement 
with the CVF data is even better. 

Fig. |fj convincingly illustrates the importance of the 
hot spots and of PAHs for the mid IR spectrum of a galac- 
tic nucleus. We mention that the contribution of the very 
small grains is almost negligible; including them hardly 
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Fig. 6. To demonstrate the significance of the hot spots and of the PAHs for the mid IR spectrum, the model of Fig. 5 
(dashed) is compared to two large-grain-only models. The first has hot spots (full line) and indicates the continuum 
underlying the PAH emission; the second is without hot spots and the dust is only heated by the ISRF. Note the 
logarithmic scale. 



changes anything. Fig. || also suggests that at 6/im, the 
wavelength where we obtained a map and performed the 
polarisation measurements, about 50% of the flux comes 
from large grains. Therefore the proposition made at the 
end of Sect. 3 that the polarisation is due to emission by 
large grains is likely. 

7. Conclusion 

ISO observations of the inner 25" x 25" nucleus of 
NGC1808 show at 6 /jm up to 20% and at 170 /jm a net po- 
larisation of 2.5%. The mid IR polarisation map is the first 
of its kind of an extragalactic source. To explain the mid 
and far IR polarisation we rule out synchrotron radiation, 
scattering from large grains, PAH or small grain emission. 
Emission by large (> 100 A), non-spherical grains, aligned 
on large scales (~ 500 pc) by uniform magnetic fields is 
proposed as the major mid and far IR polarisation driver. 
ISOCAM mid IR spectroscopy revealed a multitude of 
emission bands. We report the discovery of PAH features 
situated longward of 13^m which are detected for the first 
time in a galaxy. Based on this observation, we derive 



integrated absorption cross sections of astronomical PAH 
which can be used for future dust studies. 

The polarimetric and spectro-photometric data are 
consistently explained by a radiative transfer model which 
considers small localized regions of warm dust in the im- 
mediate vicinity of early type stars (hot spots). We demon- 
strate that this local heating gives rise to a significant 
contribution to the mid IR continuum by large grains. We 
find that the emission of small graphites or small silicates 
is negligible. 

In addition, we conclude that NGC 1808 is not a star- 
burst since in our model only 10% of the total luminosity 
comes from OB stars. This is supported by the observation 
that the optical depth in NGC 1808 is more than a factor 
5 lower than in genuine starbursts like M82 or NGC 253. 
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